Giant Spin-Hall Effect induced by Zeeman Interaction in Graphene 
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We propose a new approach to generate and detect spin currents in graphene, based on a large 
spin-HaU response arising near the neutrahty point in the presence of external magnetic field. Spin 
currents result from the imbalance of the Hall resistivity for the spin-up and spin-down carriers 
induced by Zeeman interaction, and do not involve spin-orbit interaction. Large values of the spin- 
Hall response achievable in moderate magnetic fields produced by on-chip sources, and up to room 
temperature, make the effect viable for spintronics applications. 
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The spin-Hall effect (SHE) is a transport phenomenon 
resulting from coupling of spin and charge currents: an 
electrical current induces a transverse spin current and 
vice versalHQ- The SHE offers tools for electrical ma- 
nipulation of electron spins via striking phenomena such 
as current-induced spatial segregation of opposite spins 
and accumulation of spin at the boundary of current- 
carrying sample 0, 0] ■ AH SHE mechanisms known to 
date rely on spin-orbit interaction. The two main va- 
rieties of SHE — intrinsic SHE and extrinsic SHE — 
arise due spin-orbit terms in the band HamiltonianQ and 
spin-dependent scattering on impurities respectively. 

Single layer graphene has emerged recently as an at- 
tractive material for spintronics that features long spin 
diffusion lengths 0, gate tunable spin transport 0,0], and 
high-efficiency spin injection^]. However, to realize the 
full potential of graphene, several issues must be ad- 
dressed. First, the measured spin lifetimes are orders of 
magnitude shorter than theoretical predictions 0- 11 1 call- 
ing for identifyin g an d controlling extrinsic mechanisms 
of spin scattering|10l- 14|. Second, the low intrinsic spin- 
orbit coupling values 8|, llSj render the conventional SHE 
mechanisms ineffective, depriving graphene spintronics of 
a crucial control knob for spin transport. 

Here we outline a new approach to generate and probe 
spin currents in graphene, based on a SHE response in the 
presence of magnetic field that does not rely on spin-orbit 
interaction. Spin currents are generated by the combined 
effect of spin and orbital coupling to magnetic field. The 
Zeeman splitting lifts the up/down spin degeneracy and 
imbalances the Hall resistivities of the two spin species 
(see Fig IT] inset). This leads to a net transverse spin 
current in response to an applied charge current. The re- 
sulting SHE response, called below ZSHE for brevity, is 
an essentially classical effect that offers a robust and effi- 
cient way to generate spin currents electrically in a wide 
range of temperatures and magnetic fields. The ZSHE 
response is sharply enhanced near the Dirac point (DP). 
This makes the effect viable for spintronics applications. 

The enhancement at the DP, which results from special 
transport properties of the Dirac fermions, is illustrated 
in Fig.[TJ Transport is unipolar at high doping from the 
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FIG. 1: Spin-Hall response induced by an external magnetic 
field in graphene in the absence of spin-orbit coupling. The 
SHE coefficient Ssh, Eq.Q, peaks at the Dirac point (DP). 
Spin currents at the DP originate from the the imbalance of 
Hall resistivities for spin up and down due to Zeeman split- 
ting Ez (inset, red and blue curves). Steep behavior of p^y 
leads to large imbalance in the spin-up and spin-down Hall 
response at chemical potentials < A/x. Large values Ssh 
can be reached already at moderate field strengths and high 
temperatures, Ea. (|14|l . Parameters used: B = 1 T, disorder 
broadening 7 = 100 K, electron-hole drag coefficient t) = 2.3?i. 

DP, dominated by carriers of one type, with p^y following 
the standard quasiclassical expression. 
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Transport near the DP is bipolar, which produces smear- 
ing of the 1/n singularity in p^y by the effects of two- 
particle scattering as well as disorder. This leads to a 
steep linear dependence in p^y (n) at the DP (Fig. [T] in- 
set), which is also seen in experiment (Fig.[3|). The large 
values of dpxy/dn, despite the smallness of the Zeeman 
splitting, can yield giant ZSHE response. 

The conventional SHE is described by the spin-Hall 
conductivity which relates transverse spin current and 
the electric field [H, • To identify the relevant quantity 
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FIG. 2: (a) Schematic for spin accumulation in the SHE 
regime. An electric current in a graphene strip drives trans- 
verse spin current, resulting in spin density build-up across 
the strip, Eq.(|3]). (b) Generation and detection of spin cur- 
rent in the H-geometry. Electric current passed through the 
region of local magnetic field drives spin current along the 
strip. Voltage generated via inverse SHE is detected using 
probes 3, 4. Hanle-type oscillation due to spin precession can 
be induced by external magnetic field applied in-plane. 

for ZSHE, we consider spin accumulation in a simplified 
situation when the two spin species are independent, each 
described by its own conductivity tensor. For a strip of 
width w carrying uniform current, shown in Figj^lja), the 
transverse gradients of electrochemical potential for each 
spin projection are 



^1 = (2) 
en/ pL 



with the up/down spin concentrations n-f^i^i) and the den- 
sities of states ^'f(4.). Ignoring spin relaxation, we esti- 
mate spin density at the edge = — rij^ as 
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with Ez the Zeeman splitting (for full treatment see Ap- 
pendix A). Here we used the smallness of Ez compared 
to the DP smearing Ap. (see Fig. [J) to express 9sh as a 
derivative with respect to p. Our analysis shows that 
the quantity 0sH plays a role identical to the ratio of the 
spin- Hall and ohmic conductivities ^sh = '2as}i/o'xx hr 
the conventional SHE. We will thus refer to 6'sh as the 
SHE coefficient. 

For realistic parameter values, Eq.® yields large 6'sh 
at the peak (see Fig[T|). For B = IT, using disorder 
strength estimated from mobility in graphene on a BN 
substrate, 7 « 100 K (see Eq.^), we find 6'sh = 0.1. 
This is more than two orders of magnitude greater than 
the SHE values in typical spintronics materials with spin- 
orbit SHE mechanism. Say, we estimate ^sh ~ 5 • 10^** 



from the spin and charge resistance measured in InGaAs 
system [1]. The 'giant' values 0sh are in fact to be ex- 
pected, since the ZSHE can be viewed as a classical coun- 
terpart of the SHE at /cbT < Ez discussed in Refs. 3, 12 1 
characterized by quantized ctsh = 2e^//i. 

Large Sgn values result in 'giant' spin accumulation. 
From Eq.([3|). taking 6'sh = 0.1 and the density of states at 
disorder-broadened DP ^'-|-(4,) ~ %/ An/TrSwo (with density 
inhomogeneity An sa lO^^cm"^ typical for graphene on 
BN substrate [l8[), and using £ = 10 V//zm (which corre- 
sponds to maximum current density in graphene [l9| ) , we 
estimate rig at the edges of a 2 /xm-wide graphene strip: 



n, « 3 • 10^° cm"^ 
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Such large densities, which exceed the DP width An, can 
be easily detected by spin-dependent tunneling. The esti- 
mate (|4]) is also five orders of magnitude greater than the 
spin accumulation per atomic layer observed in a three- 
dimensional GaAs [J], n^ ~ 5 • 10^ cm~^. 

Another attractive feature of the ZSHE is that it can 
enable local generation and detection of spin currents. 
Permanent micromagnets can generate fields up to IT 
concentrated to regions of size ~ 0.5 pm [i^] (fields up to 
1.4 T are achievable using widely available Neodymium 
Boron magnets). State-of-the-art microelectromagnets 
have similar characteristics [2lj. In an H-geometry, pic- 
tured in Fig. [2][b), spin currents can be generated on 
one end of graphene strip and detected on the opposite 
end. External B field, applied in-plane or at an angle to 
the graphene sheet, can be used to induce spin precession 
which will manifest itself in Hanle-type oscillations of the 
voltage measured between probes 3, 4. This setup can 
serve as an all-electric probe of spin currents 3 22, 23j . 

To model the dependence of ^sh on B, T and disor- 
der, we em ploy the quantum kinetic equation approach 
of Refs. (23. l25l|. For a spatially uniform system, we have 
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^xB)^:^ = St[/e(p),Mp)], (5) 



where /^(h) (p) is the distribution function for electrons 
and holes, and Qe — —qn = e. To describe transport near 
the DP, it is crucial to account for the contributions of 
both electrons and holes. The collision integral describes 
momentum relaxation due to two-particle collisions and 
scattering on disorder (24. 25 1. 

The approach based on Eq.® is valid in the quasi- 
classical regime, when particle mean free paths are long 
compared to wavelength. This is true when the colli- 
sion rate is small compared to typical particle energy, 
which requires weak disorder 7 <^ k^T, where 7 is de- 
fined in Ea. ([TU)) . and weak effective fine structure con- 
stant a = e^/hvQK ^ 1 (k is the dielectric constant). 

The kinetic equation (15]) can be solved analytically 
in the limit of small a |24l [25 1. Rather than pursuing 
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this route, we fohow Ref.[26[ to obtain transport coeffi- 
cients from the balance of the net momentum for differ- 
ent groups of carriers, electrons and holes, taken to be 
moving independently. We use a simple ansatz 
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, £p = vo\p\, (6) 



where fie = —fJ-h are the chemical potentials of electrons 
and holes. The quantities ag and a/j, which have the di- 
mension of velocity, are introduced to describe a current- 
carrying state. This ansatz corresponds to a uniform 
motion of the electron and hole subsystems, such that 
the collision integral for the e-e and h-h processes van- 
ishes (as follows from the explicit form of the collision 
integral given in Ref. [l^). Thus only the e-h collisions 
contribute to momentum relaxation, resulting in mutual 
drag between the e and h subsystems. 

Eq.® yields coupled equations for ensemble-averaged 
velocities and momenta of different groups of carriers ([6|) : 
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r/^rv(V,-V,0, (7) 



where i, i' label the e and h subsystems with different 
spins. The ensemble-averaged scattering times rf'^^, the 
carrier densities rii , and the electron- hole drag coefficient 
77, describing collisions between electrons and holes, are 
specified below. 

The quantities V^, are proportional to each other, 
P; = m^Vj. An explicit expression for as a func- 
tion of T, fi can be found by expanding the distribution 
functions ^ to lowest non- vanishing order in ae(;j) : 

_ 1 /rf^pp.Va./.,(p) _ 1 /rf^ppgg.,(p) 

where 5i(p) /i(p)(l - /»(?))• The integrals over p, 
evaluated numerically, give the effective mass as a func- 
tion of T and fi. At charge neutrality, setting iJie(h) = 0, 
we find ruT = ^^k^T/vl « Z.i^k^T/vl. 

The times rf^^ and carrier densities rii in ([7|) are ex- 
pressed through the distribution function ^ with a; = 0: 
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where t (e) is the transport scattering time, Eg.lfTUl). 
and the factor of two accounts for valley degeneracy. 

We pick the model for disorder scattering to account 
for the experimentally observed linear dependence of con- 
ductivity vs. doping, cr = where is the mobility 
away from the DP. This is the case for Coulomb impuri- 
ties or strong point-like defects, such as adatoms or va- 
in both cases the scattering time has an 
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approximately linear dependence on particle energy. 



kl>7 



ft|£|/7^ 7 = VQ^Jeh/fi^ (10) 



where the disorder strength parameter 7 is expressed 
through mobility. The valueu, = 6 • 10^ cm^/V • s mea- 
sured in graphene on BN [27[ yields 7 « 120 K. Similar 
values are obtained from the Paj^^-based DP width. Tak- 
ing An sa IQi" cm-2 [l3l , we find 7 ~ hvoVA^ ~ 100 K. 

To obtain pty'\ we solve Eq.©, accounting only for 
the drag between electrons and holes of the same spin. It 
can be shown that including the drag between species of 
opposite spin does not change the overall behavior of the 
transport coefficients and SHE (see Appendix C). Eq.([7]) 
can be conveniently analyzed using complex- valued quan- 
-|- iVy^ giving complex resistivity 
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Here 7^ — — ifli, with Hi — qiB/niiC the cyclotron 
frequency. 

As a sanity check, we consider the behavior at charge 
neutrality. Setting Ue ~ Uh, = rrih, etc., gives p^x 
which is a sum of the Drude-Lorentz resistivity and the 
electron- hole drag contribution analyzed in Refs. 13, 25| . 
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and p'xy' = 0. Here riT is the density of thermally ac- 
tivated electrons (holes) at the DP, having fixed spin 
projection. Disorder scattering (first term) dominates 
at low temperatures T < = l^/^Jv (at B = Q), while 
electron- hole drag (last term) dominates at T ^ T*. 

The value for the electron-hole drag coefficient 77 can 
be obtained by matching the last term in Eg. (fT2| ) . di- 
vided by 2 to account for spin, to the analytic result 
Pxx ~ SAha^/e^ [23, Ell- We evaluate a using the effec- 
tive dielectric constant k = -I- f ~ 6, which ac- 
counts for screening by substrate and for intrinsic screen- 
ing in the RPA approximation. Taking Eq « 4 for BN 
substrate [l^l, yields a ~ 0.37, giving 77 sa 2.3ft. 

The dependence of transport coefficients on T, B and 
carrier density n, predicted from Ea.(|lip. can be directly 
compared to experiment. Fig. [3] shows pxyiji) measured 
in graphene on BN, on samples similar to those described 
in Ref. [ist . The modeled Pxyiji) captures the main fea- 
tures of the data: the l/n dependence at large n and 
a steep linear dependence near the DP. The linear re- 
gion broadens with temperature at T > 7. The peak 
in Pxxin) features similar thermal broadening (see Ap- 
pendix C). The SHE coefficient, found from Eq.([2]), is 
plotted in Fig[TJ 

We now explore the behavior of transport coefficients 
near the DP, making estimates separately for T ^ and 
T ^ T*. This can be conveniently done using an interpo- 
lation formula Tj^'^(p,T) = mi{p,T)vQh/j^ which links 
the ensemble-averaged scattering time ^ and the effec- 
tive mass ([S]) in the entire range of T and fi of interest. 
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Since our SHE mechanism does not rely on the rela- 
tivistic dispersion of excitations, it can also be realized in 
other zero-gap semiconductors, in particular graphene bi- 
layer. It also applies, with suitable modifications, to the 
valley degrees of freedom in graphene. It was predicted 
that a (non-quantizing) magnetic field can lift valley de- 
generacy and produce a Zeeman-like valley splitting [s^ . 
This will imbalance the Hall resistivities and result in a 
valley-Hall effect of a magnitude similar to the SHE. 

We thank D. Goldhaber-Gordon, L. M. K. Vander- 
sypen and M. Soljacic for useful discussions and acknowl- 
edge support from Naval Research Grant NOOO 14-09-1- 
0724 (LL). 



FIG. 3: Measured pxy{n) for a high-mobility graphene sample 
on BN substrate at T = 250 K. The dependence follows the 
quasiclassical formula ([l]) away from the DP, and is linear with 
a steep slope near the DP. Inset: Results for Pxyin) obtained 
from the two-carrier model, Eas.(f7)l. (|ll[l . for disorder strength 
7 = 180 K found by fitting the min/max distance in measured 
Pxy for B = 1 T. Other parameters: r? = 2.3ft, T = 250 K. 

We find the slope of p^y at the DP by expanding 
Eq.(|lip in small n = Ue — Uh (see Appendix D for full 
treatment). The result, which simplifies in each of the 
regimes T ^T^, and T ^T^., can be described by a single 
interpolation formula as 
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where only terms first-order in B have been retained. 

The SHE coefficient, Eq.Q, found by combining the 
results ([T^ and ([T^. and using thermally broadened 
density of states at the DP derived in Appendix B, 
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where Eq is the cyclotron energy. The functional form 
is the same in both regimes, dsn cx B^ /T, with different 
prefactors \t>t, = 241n2/7r^ and \t<t, = 121n2/7r^. 
The 1/T growth of Osn saturates at k^T k, 7, reaching 
maximum value 0sH,max ~ h^E^Ez/^'^ 



We expect suspended graphene [28|, l29j to feature an 
even stronger SHE than graphene on BN. Using typical 
mobility = 2 • lO^cm^/Vs [i^, we estimate 7 ~ 65K, 
whereas the temperature dependence of the conductivity 
at the DP [H] yields 7 10 K. For either value of 7, 
Eq.dTH) predicts very large values 9sh at the DP. 

Based on these estimates, we expect strong SHE re- 
sponse already at moderate fields i? ^ 1 T. Besides spin 
accumulation and locally tunable SHE response, which 
was discussed above, SHE can also manifest itself in a 
non-zero Hall voltage in response to spin-polarized cur- 
rents injected from magnetic contacts. 
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APPENDIX A: SPIN ACCUMULATION 

Here we analyze spin accumulation in the ZSHE 
regime. Spin accumulation in two-dimensional electron 
gasses, arising due to weak SHE induced by spin-orbit 
interaction, was considered in several papers, in particu- 
lar: S. Zhang, Phys. Rev. Lett. 85, 393 (2000). There 
are several new aspects in our problem that warrant spe- 
cial treatment. Specifically: (i) In the presence of exter- 
nal magnetic field charge transport is described by p^x 
and pxy separetly for each spin projection; (ii) Unlike 
the semiconductor case, the expected SHE response is 
not necessarily small, ^gu ^ 1. 

We consider spin accumulation in the strip geome- 
try, —w/2 < y < w/2, with current driven along the 
strip, such that external electric field £ is applied along 
X direction. To solve the transport problem, we intro- 
duce electrochemical potentials for each spin projection, 
'^t(i) = + ^T(^)/'^t(-l)' where (j) is the electric poten- 
tial, n^[i) is the deviation of the local density for spin up 
(down) from its equilibrium value, and i^^d) is the den- 
sity of states for the two spin projections. The equations 
for the current density are given by. 



For a narrow strip, w <^ is-, spin relaxation can be ig- 
nored. In this case, the expression for ns{y) agrees with 
the estimate ([3]) in the main text. 

To estimate the numerical value of the spin density 
at the DP, where it is maximum, we make several as- 
sumptions. First, we assume that the density of states 
is disorder-broadened at the Dirac point, i? = v-^ = !v 
SI i , where 7 is the disorder strength parameter defined 



j|(r) = (7'^V^t(r), j^(r) = a^W(pi{r 



(15) 



where f and J, label carriers with up-spin and down-spin. 
Here (t^ and &^ are 2x2 matrices describing the longi- 
tudinal and Hall conductivity. The continuity equation 
for each spin projection can be written as 

Vjt(r) = --fs{n^{r)-ni{r)), Vjii^) = -7,(n^(r)-nt(r)), 

(16) 

where 7^ is the rate of spin relaxation. The terms with 
time derivative are omitted, as appropriate for the DC 
transport. 

The transport equations should be supplemented by 
the condition of electro- neutrality, n-|-(r) = —ni{r). We 
solve transport equations in a strip geometry, where cur- 
rent is driven by an electric field parallel to the strip, 
with the boundary conditions of zero normal current at 
the strip edges, jy{y = ±|w) = 0. We obtain the spin 
density rig = n-f — profile across the strip 
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where v = 2/ [v^^ + v'^^) is the average density of states 
for the two spin projections, £ is the electric field, is is 
the spin relaxation length 



in Eq. ([TO| of the main text. For an estimate we use the 
value 7 ~ 100 K, which corresponds to density inhomo- 
gencity An w 2v{Q)^ sa 10^°cm~^ typical for graphene 
on BN substrate. Second, we will assume w ~ 2fim, and 
for the SHE coefficient, we will assume ^sh = 0.1, as es- 
timated in the main text. Taking £ = 10V//im, as in 
Ref. [1^, for the spin density at the edge we find 



Usiy = ±w/2) « 3 • 10 



10 cm-2. 



(19) 



This is at least five orders of magnitude the spin density 
per atomic layer observed in GaAs Q , ~ 5 x 10^ cm~^. 
Such giant spin accumulation is due to the larger value 
of 0SH in graphene, and due to the fact that graphene 
can sustain large current densities and electric fields. 



APPENDIX B: THE CARRIER DENSITY AND 
THE DENSITY OF STATES 

We start with deriving a closed form expression for 
the carrier density, evaluated separately for the elec- 
trons and holes with a fixed spin projection: ni(/^,T) = 
2 / (27r^a /i(P; T), where the distribution function fi is 
defined by Eq.® of the main text with = 0. Integrat- 
ing over p, we find the dependence of ^g^/j) on tempera- 
ture and chemical potential as power series in g = e'^/^'^'^ 
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Although the series converge only for — 1 < g < 1, the 
final expression is valid for both 5 < 1 and g > 1. At the 
DP, setting /i = 0, we obtain 



^e(h)(/^ = 0) = 
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This quantity, denoted as ht in Eq. ()12|) of the main text 
and elsewhere, gives the density of thermally excited car- 
riers (electrons or holes) with a fixed spin projection. 

Next, we analyze the density of states f = dn/dp. The 
density of states per one spin projection is represented by 
a sum of an electron and a hole contribution 
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where the minus sign reflects the fact that electrons and 
holes have opposite chemical potentials, /ie = ~^J'h = 
We use the above expression for ne(/j) to find 

= (ln(l+5)+ln(l+.g-i))-M_ (23) 
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This expression describes temperature-broadened density 
of states: 
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with a crossover temperature k^T = |/.i|/(21n2). 



(25) 



APPENDIX C: THE EFFECT OF DRAG 
BETWEEN CARRIERS OF OPPOSITE SPIN 

Our analysis of transport coefficients in the main text 
was based on a simplified model which neglected drag be- 
tween carriers of opposite spin. To evaluate the accuracy 
of this approach, here we discuss a more general model 
which accounts for drag between carriers of either spin. 
To emphasize this difference, we will use notation fj for 
the drag coefficient instead of rj used in the main text. 
As we will see, the more general approach predicts an es- 
sentially identical behavior of transport coefficients, and 
a qualitatively similar behavior of the SHE coefficient. 

We consider coupled transport of four carrier species. 
For simplicity, we will take the drag coefficient values to 
be the same for all carrier species, a reasonable approx- 
imation at high kBT. Ensemble-averaged velocities and 
momenta in the presence of electric and magnetic fields 
are governed by Eq.Q of the main text, which we dupli- 
cate here for reader's convenience. 



,(V,-V,0. (26) 



The quantities rii, rrii, are defined as in the main 
text, with the chemical potential /x replaced by 
for spin up (down) electrons. The value of the dra g co - 
efficient is found by matching analytic results[2^ l25j. 
Because Eas. ipS]) account for drag between all carrier 
species, not just those with parallel spins, the value fj 
found below is different from the one obtained in the 
main text. 

The transport coefficients are determined as follows. 
First, for given values fj, and T, we obtain rii, uii by 
numerically evaluating integrals in the relation 



r/^(P,/^,T), 



(27) 
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FIG. 4: The longitudinal resistivity, obtained from the model 
(|26p which takes into account drag between carriers of either 
polarity and spin, shown for several temperatures. Param- 
eters used: B = OT, 7 = lOOK, r? = 1.15ft. The peak in 
Pxx gets broadened as the temperature increases, as pxx{0) 
approaches the high-temperature limiting value, Eq.((3l]). 



and in 



Eq.i 



For ensemble-averaged scattering time 



use the interpolation formula T^(fi,T) 



■mi{fi,T)vQh/j'^ which links the ensemble-averaged scat- 
tering time ([9]) and the effective mass ^ in the entire 
range of T and n of interest. 

After that, we solve the four coupled linear equations 
(|26p and find the currents of the spin-up and spin-down 
electrons, 

h = "-eteVet - Hh^eYh^, = UeieVei - nhieVhi- 

(28) 

Then the charge and spin currents are expressed through 
jf and jj. as follows. 



Jc— Jt+Jj-j Js— Jt~J4,J 



(29) 



These expressions can be used to calculate the charge and 
spin conductivity tensors. We obtain the dimensionless 
SHE coefficient by evaluating the ratio of the transverse 
spin current and the longitudinal charge current. 



(30) 



where the transverse and longitudinal components are 
taken with respect to the electric field E. 

For charge transport, this model leads to the behav- 
ior of transport coefficients which is essentially identical 
to that obtained from the simplified model used in the 
main text, Eq. ([TT]) . albeit with a doubled value of the 
drag coefficient, 77 2fj. To see this, we note that the 
forces qe(h) + 'j'^ ^ ^) that drive transport are the 
same for the up-spin and down-spin carriers. Further, 
the quantities Ui, im, rf"* coincide for spin up and down 
when Ez = 0. Therefore, in the limit — > 0, the rela- 
tions (1^51) can be satisfied by v| = vt, vt = v^. Tak- 



ing this into account and eliminating variables for one 
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FIG. 5: The SHE coefficient ^sh, Eq.((30l), obtained from tlie 
model (|26p which takes into account drag between carriers of 
either polarity with parallel an opposite spins. Parameters 
used: B = 1 T, 7 = 100 K, 77 = 1.15ft. Similarly to ^sh, illus- 
trated in Fig. [1] ^SH exhibits a sharp peak at the DP, which 
gets more pronounced as T is lowered. The quantities ^sh 
and 6'sH are identical at low T, when transport is dominated 
by scattering on disorder. At high T, when transport is dom- 
inated by the effects of electron-hole drag, the peak in ^sh is 
significantly broader than the peak in 6'sh. 



spin projection wc obtain equations for the other pro- 
jection which are identical to the equations in the main 
text up to a replacement ry — >■ 2ry. Small but finite Ez 
changes the transport coefficients, however the difference 
between the spin-up and spin-down remains small as long 
as Ez < fcBr,7. 

Following the same reasoning as in the main text, we 
fix the drag coefficient value at the half the drag coeffi- 
cient used in the simplified model: 2fi = 2.3/1. We obtain 
the density dependence Pxy{n) which is very similar to 
that shown in Fig. [3] of the main text. Behavior of the 
longitudinal resistivity Pxxin) at i? = is illustrated in 
Fig. m The peak in p^x at the DP becomes lower and 
broader as the temperature is increased. The resistiv- 
ity at the DP decreases, at high temperatures /jbT ^ 7 
saturating at the value 



1 



Px 



,t(i) 



5m, 



(31) 



where the factor 1/2 is introduced to convert the resis- 
tivity for one spin projection, given by Eq.(fT2|) of the 
main text, to net resistivity. The results in Fig. |4] show 
that Pxx{0) is very close to this limiting value already at 
T ~ 300 K, indicating that transport at these tempera- 
tures is dominated by electron- hole collisions. 

The behavior of the SHE coefficient ^sh, Eg.lpO]). 
shown in FiglU is overall similar although not identical 
to the behavior of ^sh found in the main text. Near the 
DP, ^sH reaches large values, similar to those of Osn (see 
Fig. [T] in the main text). The peak in ^sh is, however. 



significantly broader. Furthermore, ^SH exhibits a mono- 
tonic decay away from the DP; in contrast, Ogn exhibits 
a sign change. We believe the sign change to be pecu- 
liar for the two-component model which ignores the drag 
between spin-up and spin-down carriers. 

The SHE coefficient £^sh determines physical observ- 
ables, such as spin accumulation density, which is given 
by Eq.Q of the main text with a substitution 9sii — >■ ■Csh- 
The two models therefore predict the same values of Us 
at the DP. 



APPENDIX D: ANALYTIC ESTIMATES FOR 
THE SHE COEFFICIENT 

Here we derive the expression for dpxy/dn, Ea. (|T3l) . 
and an estimate for the SHE coefficient, Eq. ([T4| of the 
main text. First, we consider the high-temperature 
regime, T 3> T*, where transport coefficients are dom- 
inated by the effects of electron-hole drag. In this limit, 
the first term in the numerator of Ea. ((TT|) . which is small 
compared to the other two terms, can be neglected. Fur- 
thermore, the last term in the denominator, which is 
quadratic in density, cannot affect the slope of pxy at 
the DP, and we can drop it as well. Then Eg. (ITT]) is 
simplified as 



[Pxx + ipxy) 



t(i) 



-le + ^7/. 



52 ^7h + ^7e 



(32) 



Using the relations 7^ = — iVli and the interpolation 
formula rf'^ = rriiv'^h/^'^ , after simple algebra we obtain 

(Pxx-tipxy) 2~27 1 \ / 2fc , ~15 ]~ ' y'^'^) 

7^(rie 4- nhjlv^n + leBn/c 

where n = — is the carrier density. Then, using 
the fact that Tig + = 2nT + 0{n^), and expanding to 
linear order in B, we obtain 



g2 ' P'^y 



h^vl B 
T2 nxec 



(34) 



The second relation gives the slope of the Hall resistivity 
at the DP, Eq.(Il3l) of the main text (case T >T^). 

In the low-temperature regime, T <^T^, we can neglect 
the second and third terms in the numerator of Eq. pT|) ,as 
well as the third term in the denominator, which gives 



{Pxx ~t~ ipxy ) 



t(4) 



lelh 



^7e 

rah ' 



(35) 



ifi, and the 



Once again using the relations 7^ = 

interpolation formula rf^^ = m.iv'^h/^'^ , we rewrite the 
above expression as 



{Pxx ~t~ ipxy) 



1 



"f^ijie + nh)/hvQ + ieBn/c 



(36) 
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where we have neglected the term proportional to fle^h 
in the numerator, which is of the order B^. Expanding 
to linear order in ri, wc obtain 

_ 7' _ _B n_ 

Using the relation nx = ('fef") ' obtain the slope 



of pxy at the DP, Eq.dlS]) of the main text (case T < T*). 

Finally, we obtain the SHE coefficient at the DP, de- 
fined in Eq.(I3]) of the main text. This is done by combin- 
ing the results for dpxy/dn, which we have just derived, 
with the result for thermally broadened density of states 
at the DP, Eg.pS]). and the above results for p^x at the 
DP. This gives Ea. ([ni) of the main text. 



